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The structure of the concentration field of a decaying 
substance produced by chemical sources and advected by a 
smooth incompressible two-dimensional flow is investigated. 
We focus our attention on the non-uniformities of the Holder 
exponent of the resulting filamental chemical field. They ap- 
pear most evidently in the case of open flows where irregular- 
ities of the field exhibit strong spatial intermittency as they 
are restricted to a fractal manifold. Non-uniformities of the 
Holder exponent of the chemical field in closed flows appears 
as a consequence of the non-uniform stretching of the fiuid 
elements. We study how this affects the scaling exponents 
of the structure functions, displaying anomalous scaling, and 
relate the scaling exponents to the distribution of finite-time 
Lyapunov exponents of the advection dynamics. Theoretical 
predictions are compared with numerical experiments. 

PACS: 47.52.-fj, 05.45.-a, 47.70.Fw, 47.53.-fn 



I. INTRODUCTION 

Mixing in fluids plays an important role in nature and 
technology with implications in areas ranging from geo- 
physics to chemical engineering The phenomenon of 
chaotic advection - intensively investigated during the 
last decade - provides a basic mechanism for mixing 
in laminar flows ||^. Briefly stated, chaotic advection 
refers to the situation in which fluid elements in a non- 
turbulent flow follow chaotic trajectories. Advection by 
simple time-dependent two-dimensional flows falls gener- 
ically under this category. Stirring by chaotic motion, 
with its characteristic stretching and folding of material 
elements, is able to bring distant parts of the fluid into 
intimate contact and thus greatly enhances mixing by 
molecular diffusion acting at small scales. 

Mixing efficiency becomes specially important when 
the substances advected by the flow are not inert but 
have some kind of activity. By 'activity' we mean that 
some time-evolution is occurring to the concentrations 
inside advected fluid elements (produced by chemical re- 
actions, for example). For definiteness we will use terms 
such as chemical fields and chemical reactions, but biolog- 
ical processes, occurring for example when the advected 
substance is living plankton, can be described formally 
in the same way. The interaction between the stirring 
process and the chemical activity can result in complex 



patterns for the spatial distribution of the chemical fields, 
which in turn greatly affect the chemical processes [^1). 
In addition to the impact on its own chemical dynamics, 
the spatial inhomogeneities may have important effect 
on other dynamical processes occurring in the fluid (for 
example in the behavior of predators seeking for the ad- 
vected plankton |^). An understanding of the structure 
of these spatial patterns is thus valuable. 

Previous theoretical work concentrated on the tempo- 
ral evolution of the total amount of chemical products in 
specific reactions such a.s A + B 2C or A + B ^ 2B 
||. In the same type of reactions were studied in 
open flows. In a previous paper Q some of us con- 
sidered a class of chemical dynamics characterized by 
a negative Lyapunov exponent in the presence of (non- 
homogeneous) chemical sources. Under such chemical 
processes, reactant concentrations present a tendency to 
relax towards a local-equilibrium concentration (the fixed 
point of the local chemical dynamics). This tendency is 
disrupted by the advection process, which forces fluid 
elements to visit places with different local-equilibrium 
states. Depending on the relative strength of chaotic ad- 
vection and relaxation the resulting concentration distri- 
bution can be smooth (differentiable) or exhibits char- 
acteristic filamental patterns that are nowhere differen- 
tiable except in the direction of filaments aligned with 
the unstable foliation induced in the fluid by the chaotic 
dynamics. The mechanism for the appearance of these 
singular fllaments is similar to the one producing singu- 
lar invariant measures in dynamical systems , although 
here it is affected by the presence of the chemical dy- 
namics: stretching by the flow homogeneizes the pattern 
along unstable directions, whereas small-scale variance, 
cascading down from larger scales, accumulates along the 
stable directions, producing diverging gradients. 

The strength of the singularities of the concentration 
fleld C (r) can be characterized by a Holder exponent a 

\SC{vo:.Sr)\ = |C7(ro + 6r) - C(ro)| ~ |<5rr, \Sr\ ^ . 

(1.1) 

If the field is smooth (differentiable) at ro, a = 1, 
while for an irregular rough (e.g. filamental) structure 
< a < 1. In [|| we focused on the existence of a 
smooth-filamental transition as time-scales of the system 
are varied, and also obtained the most probable (bulk) 
value of the Holder exponent. Note, however, that the 



1 



Holder exponent defined by (1.1) is a local characteristic 
of the field, whose value may depend on the position Tq . 
In this Paper we concentrate on such non-uniformities of 
the filamental chemical field and study how this affects 
scaling properties of quantities involving spatial averages, 
which are the more convenient quantities to be observed 
in experiments. 

In Section |l| we review the results presented in 
namely the smooth-filamental transition and the domi- 
nant value of the Holder exponent in closed flows. Then 
we consider the sa me problem for the case of mixing by 
open flows (Sect. ^I^ . In this case the necessity of a 
multifractal description becomes manifest, and this mo- 
tivates the development of a quantitative characteriza- 
tion of the filamental structures in terms of structure 
functions. This is presented in Sect. Scaling expo- 
nents appear to be related to the distribution of finite- 
time Lyapunov exponents. We conclude the Paper with 
a summary and discussion. 



II. THE SMOOTH-FILAMENTAL TRANSITION 

We consider the flow as externally prescribed, thus ne- 
glecting any back influence of the chemical dynamics into 
the hydrodynamics (the advected substances are chemi- 
cally active but hydrodynamically passive). In this con- 
text, the general continuum description of chemical reac- 
tions in hydrodynamic flows is given by sets of reaction- 
advection-diffusion equations. They involve in general 
multiple components and nonlinear reaction terms. Ref- 
erence ^ considered the situation in which the chemical 
kinetics is stable, i.e. there is a local-equilibrium state 
at each spatial position, determined by the sources and 
the reaction terms, so that concentrations of fluid par- 
ticles visiting that position tend to relax to the local- 
equilibrium value. Mathematically this corresponds to 
the negativity of the Lyapunov exponents associated to 
the chemical dynamical subsystem. It was shown in H] 
that arbitrary chemical dynamics in this class can be sub- 
stituted by linear relaxation towards local equilibrium at 
a rate given by the largest (least negative) chemical Lya- 
punov exponent. Within this restriction, the multiplicity 
of components is not essential since, except for special 
types of coupling, linearization leads to simple relations 
between the different fields. 

Because of the above remarks, and with the aim of 
keeping the mathematics as simple as possible, we will 
restrict our considerations in this Paper to the simplest 
chemical evolution: linear decay, at a rate 6, of a sin- 
gle advected substance. A space-dependent source of the 
substance will also be included, to maintain a non-trivial 
concentration field at long times. This chemical dynam- 
ics can be considered either as an approximation to more 
complex chemical or biological evolutions, with maxi- 
mum chemical Lyapunov exponent —b, or as a description 
of simple specific processes such as spontaneous decom- 



position of unstable radicals, decay of a radioactive sub- 
stance, or relaxation of sea-surface temperature towards 
atmospheric values 1 10 1 . The validity of our ideas for non- 
linear multicomponcnt situations has been checked for a 
plankton model in ||ll| . 

The concentration field C(r,t), when advected by a 
incompressible velocity field v(r, t) is governed by the 
equation 



dC 
'dt 



v(r, t)-VC^ S{t) -bC + kW^C, 



(2.1) 



where k is the diffusion coefficient, b is the decay rate in- 
troduced above and S{r) is the concentration input from 
chemical sources (negative values representing sinks). We 
restrict our study to the case in which the incompress- 
ible velocity field is two-dimensional, smooth, and non- 
turbulent. Chaotic advection is obtained generically if 
a simple time-dependence, for example periodic, is in- 
cluded in v(r,t). We assume that diffusion is weak and 
transport is dominated by advection. Thus one expects 
that the distribution on scales larger than a certain dif- 
fusive scale is not affected by diffusion. Therefore we 
consider the limiting non-diffusive case k — 0. In this 
limit the above problem can be described in a Lagrangian 
picture by an ensemble of ordinary differential equations 



dr 

dt 



'dt 



Sim] - bC, 



(2.2) 



(2.3) 



where the solution of the first equation gives the trajec- 
tory of a fluid parcel, r{t), while the second one describes 
the Lagrangian chemical dynamics in this fluid element: 
Cit) = C[r^r{t),t]. 

To obtain the value of the chemical field at a selected 
point r at time i one needs to know the previous history 
of this fluid element, that is the trajectory r{t) {0 < t < t) 
with the property r{t) = r. This can be obtained by the 
integration of ( p^ ) backwards in time. Once f(t) has 
been obtained, the solution of (2.3) is 



C(r,i) = C[f(0),0]e-''*+ / S[r{t)K''^'''^ dt. (2.4) 

Jo 

One can obtain the difference at time t of the values 
of the chemical fleld at two different points r and r + 5r 
separated by a small distance Sr in terms of the difference 
C'[f(i)-f<5r(t),i]-C[f(i),<] = SC[r{t),tiSr{t)] for < t < 
t, namely: 



-fct 



SC{f, t; Sr) = (5C[f (0), 0; <5r(0)]e 

+ / (55[f(t);(5r(i)]e~''(*-*)dt 



(2.5) 



where Sr{t) {0 < t < t) is the time-dependent distance 
between the two trajectories that end at r and r -|- ifr at 
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time t, and 6S, in analogy with SC, is the difference of 
the source term at points f (t) and r{t) + Sr{t). 

Thus we have expressed the behavior of an Eulerian 
quantity SC{f^t;Sr) in terms of Lagrangian quantities, 
in particular of Sr{t). Further analysis of Eq. ( p.5| ) re- 
quires specification of the behavior of Sr. The signature 
of chaotic advection is the exponentially growing behav- 
ior of this quantity at long times. More precise state- 
ments need additional assumptions on the character of 
the flow. The simplest framework is obtained if we re- 
strict our attention to initial conditions f (0) in an invari- 
ant hyperbolic set ^ . In this case one can identify at each 
point two directions, the one in which the flow is con- 
tracting c(r) and the expanding direction e(r). They de- 
pend continuously on position (for time-dependent flows 
there is an additional explicit time-dependence that we 
do not write down to simplify the notation). Since c 
and e form at each point a vector basis which is not or- 
thonormal, it is convenient to introduce also the dual ba- 
sis (c^,e^) at each point. Chaotic advection manifests in 
which for most initial separations Sr{0) the long-time be- 
havior of Sr is \Sr{t)\ ~ e[f(0)]'f • Sr{0) e^* (for t > 0), 
where A is the positive Lyapunov exponent of the flow 
and e[f (0)]''' • Sr{0) gives the component of the initial sep- 
aration along the expanding direction. At long times, 
the direction of Sr{t) tends to become aligned with the 
expanding direction of the flow at f (i), e[f (i)]. However, 
if the initial separation is aligned with the contracting 
direction at the initial point, A should be substituted by 
A', the contractive Lyapunov exponent, and e'l' by c^. For 
incompressible flows o ne h as A' = —A. 

In order to analyze (2^) one has to consider the back- 
wards evolution. In this case typical solutions behave. 



for t < and large, as |i5r(t)| 
c[f (0)]t • <5r(0) 



c[f (0)]t • ,5r(0) 



SO that, also in this backwards dy- 
namics, close initial conditions diverge, and the difference 
will tend to become aligned with the most expanding di- 
rection of the backwards flow (the contracting of the for- 
ward flow, c[f(t)]). Again, there is a particular direction 
for the orientation of the initial condition (the contract- 
ing one in the backwards flow which is the expanding 
one in the forward dynamics) for which —A should be re- 



placed_by A. In (|2.5| ), Sr is obtained backwards starting 
from Sr a.t t = t. In this case: 



Sr{t) w c(r)t • Sr e^^'''>c[r{t)] , ts < t < t . (2.6) 



The exponential separation ( p.6[ ) holds only while the 
distance Sr{t) is not too large, and saturates when ap- 
proaching the size of the system or some characteristic 
coherence length of the velocity field. The time at which 
this happens defines t^, a saturation time. We assume 
that both the velocity field and the source S{r) have only 
large-scale structures such that their corresponding co- 
herence lengths are comparable to the system size, that 



we take as the unit of length-scales. Thus, the saturation 
time is given by tg = t + X~^\n\Sr\. 

For small Sr, Eq. (|2.5|) can be written as 



SC{f, t; Sr) w (SC[f (0), 0; 5r{Q)]e 



SS[r{t)-5r{t)]e-^^^-'Ut 
+ / 5r{t) ■ V5[f (t)]e-''(*-*)dt 



(2.7) 



(if tg > 0). We will not need to specify the dependence of 
SS on Sr for times previous tots, as long as SS remains 
bounded. Substitution of ( ^.6D in the second integral 
leads to 



SC{f, t; Sr) w (5C[f (0), 0; 5r{Q)]e 



SS[r{t);Sr{t)]e-^^^-'Ut 
+5r ■ c(r)^ / c[r{t)] ■ V S[r{t)]e^^-'''^^^-'Ut (2., 



Taking the limit (5r — > (for a finite t) leads to < 
0. Thus the first integral disappears and the first term 
can be linearized. By writing Sr = r\.\Sr\ so that n is a 
unit vector, one finds the directional derivative along the 
direction of n as 



n • VC(r, t) « n • c(r)Tc[f (0)] • VC[f (0), 0] 



=(A-b)t 



+n ■ c(r)t / c[f (t)] • V5[f (i)]e(^-^'(*-*)dt . (2.9) 



If A < 5 this derivative remains finite in the t — > oo 
limit and the asymptotic field Coo(r) = C(r,t oo) is 
smooth (differentiable) . Otherwise the derivatives of C 
diverge as ~ e'-^"''^* leading to a nowhere-differentiable 
irregular asymptotic field. The exception again is the 
expanding direction of the forward flow: when Sr points 
along that direction one should substitute in Eq. ( |2.9| ) A 
and c by —A and e. This directional derivative is always 
finite. Thus, there is at each point a direction along 
which Coo is smooth. By changing the values of b and A 
one will encounter, at 6 = A a morphological transition 
between a smooth pattern and a filamental (i.e., singular 
in all but one directions) one. It should be noted that, 
because the explicit time-dependence of the vectors c and 
e referred to before, the limiting distribution Coo will not 
be a steady field, but one following the time dependence 
of the stable and unstable directions. For time-periodic 
flows v(r,f). Coo will also be time periodic. Its singular 
characteristics however do not change in time. 

In order to characterize the singular asymptotic field 
we take the limit f — > oo for fixed finite (5r in (2. 



(5Coo(r;(fr)w / SS[r{x);5r{x)]x-^-^dx 

J|(5r|-1/A 

\Sr\n ■ c{fY J c[r{x)] ■ VS[r(x)]x^-''-'^dx (2.10) 
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where we used the change of variables e*~* ^ x. If & > A 
one finds for the dominant term in the \5v\ -^0 hmit the 
simple scaling 5C oo ^ |^r|, but when b < X: 



SCoo{r;Sr)^\Sr\ 



(2.11) 



According to (1.1) the value of the Holder exponent is 



a = min ^ ^ 1 1 



(2.12) 



This implies that if 6 < A the asymptotic chemical field 
becomes an irregular fractal object (again there is an 
orientation of Sr along which a = 1). Consequently, the 
graph of the field along a one-dimensional cut or transect, 
or the contours of constant concentration are also frac- 
tals, as they are two-dimensional sections of the whole 
Coo(r) embedded in a three-dimensional space. More 
precisely, the one-dimensional transect (along the direc- 
tion x) of the field is a self-affine function with its graph 
embedded in an inherently anisotropic space (C — x) with 
axis representing different physical quantities. Contours 
of constant concentrations, however, are self-similar frac- 
tal sets of the two-dimensional physical space {x — y). 

In Figs. |l| and H we present sna psho ts of the asymptotic 
field Coo evolving according to (2.2)-(pT^). For the flow 
we take a simple time-periodic velocity field defined in 
the unit square with periodic boundary conditions by 

Vx{x,y,t) = - t mod Tj cos(27ry) 



2C/. 



T 



Vy{x,y,t) = -^^[t mod T- — j cos(27rx) (2.13) 

where Q{x) is the Heavyside step function. In our sim- 
ulations U = 1.2, which produces a flow with a single 
connected chaotic region. The value of the numerically 
obtained Lyapunov exponent is A ~ 2.67/r. 

Backward trajectories with initial coordinates on a 
rectangular grid were calculated and used to obtain the 
chemical field at each point by using ( |2.4| ) forward in time 
with the source term S{x,y) = 1-1- sin (27rx) sin (27rj/). 
The values of the parameters used in Fig. |l| are T = 1.0 
and b = 4.0, for which the Lyapunov exponent is A = 
2.67 < 6. A smooth pattern is seen, in agreement with 
our theoretical arguments. In Fig. ^ the parameters are 
T = 1.0 and b = 0.1, so that A > 6 and a filamental 
pattern is obtained. 

The smooth-fractal transition also appears in the time- 
dependence of the concentration measured at a fixed 
point in space. This can be shown by a similar analy- 
sis for the difference 



6C{r, t; 6t) = C(r, t + 6t) - C(r, t) 



(2.14) 



instead of the spatial difference discussed above. If 
6 < A the signal C(r, t) becomes non-differentiable in 
time and can be characterized by the same Holder expo- 
nent b/X. The fact that scaling properties of the tem- 
poral signal and that of the spatial structure are the 



same -analogously to the so called 'Taylor hypothesis' in 
turbulence- can be exploited in experiments or in analy- 
sis of geophysical data. 

We conclude with some comments on the range of va- 



lidity ofEq. (2.£). The Lagrangian description (2.2 )-( 2.2) 
in which our approach is based is valid only for scales at 
which diffusion is negligible. Thus there is a minimum 
admissible value Idiffin) ~ ^/k of i5r and our calcula- 
tion should be understood as giving the gradients only 
up to this scale, fractality being washed out at smaller 
scales by the presence of diffusion. Nevertheless we think 
that, if diffusion is weak, the fractal-filamcntal transition 
will be seen at scales larger than this diffus ive s cale. For 
fixed Sr larger than the diffusion length, (2.8) remains 
valid until a time i ^ ts{Sr) that means that the diver- 
gence of the gradients will also saturate at a finite value 

Another limitation to the validity of our equations 
arises from the fact that, for most chaotic flows of phys- 
ical relevance, not all the points visited by the fluid par- 
ticle will be hyperbolic. The stable and unstable direc- 
tions in the previous discussion beco me tangent at some 
points and equations such as 
there 



(2., 



become undeflned 



More importantly, KAM tori will be present in the 
system, so that for values of r lying on KAM trajectories 
the value of the Lyapunov exponent appearing in (2.6) 
will not be positive but zero. These facts imply that our 
previous equations are not valid in a global sense. The 
main effect of the KAM tori will be to partition space into 
ergodic regions. Each region will be characterized by a 
different value of A, and there is the possibility of observ- 
ing different morphologies in the different regions. For 
the values of parameters used in Figs. |l| and ^ KAM tori 
occupy a very small and practically unobservable por- 
tion of space, so that the filamental pattern appears to 
be well described by the same Holder exponent nearly ev- 
erywhere. However, an example is given in H] in which 
filamental and smooth regions coexist separated by KAM 
tori. 

It is known that, even inside ergodic regions, the Lya- 
punov exponent is a non-smooth function of space [[l2| : 
it can differ from the most probable value A in sets of 
zero Lebesgue measure. This stresses again that our 
arguments above have a local nature. Equations such 
as ( 2.12| ) give the value of the 'bulk' or 'most probable' 
Holder exponent, that is, the one which characterizes al- 
most all points in the fluid. More generally one should de- 
fine local Holder exponents a(r) which depend on space. 
Deviations from the most probable or bulk value occur 
only in portions of the ergodic regions having vanishing 
measure. They can however have effects on the physi- 
cally observable average quantities. These effects appear 
in a specially clear manner when discussing open flows, 
the subject of the next Section. 
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III. OPEN FLOWS 

Let us now consider again the problem (^.2|)-( pT^ ) with 
a velocity field corresponding to an open flow whose time- 
dependence is restricted to a finite region (mixing region) , 
with asymptotically steady inflow and outflow regions. 
A prototype of this flow structure is a stream passing 
around a cylindrical body. If the inflow velocity is high 
enough vortices formed in the wake of the cylinder make 
the flow time-dependent in this region, while the flow re- 
mains steady in front of the cylinder or in the far down- 
stream region. We assume again that the flow is non- 
turbulent, so that the velocity field is spatially-smooth. 

Passive advection in such open flows was found to be 
a nice example of chaotic scattering 1 13 ij] . Advected 
particles (or fluid elements) enter the unsteady region, 
undergo transient chaotic advection, and finally escape 
and move away downstream on simple orbits. The time 
spent in the mixing region, however, depends strongly on 
the initial coordinates, with singularities on a fractal set 
corresponding to particles trapped forever in the mixing 
region. This is due to the existence of a non-attracting 
chaotic set (although of zero measure) formed by an infi- 
nite number of bounded hyperbolic orbits in the mixing 
region. The stable manifold of this chaotic set (or chaotic 
saddle) contains orbits coming from the inflow region but 
never escaping from the mixing zone. These points cor- 
respond to the singularities of the residence time. If a 
droplet of dye is injected into the mixing region, most 
of it will be advected downstream in a short time. But 
part of the dye will remain close to the chaotic saddle 
for very long times, and continuously ejected along its 
unstable manifold. In this way the dye traces out the 
unstable manifold of the chaotic saddle, resulting in frac- 
tal patterns characteristic to open flows |l^Jl^. Perma- 
nent chaotic advection is restricted to a fractal set of zero 
Lebesgue measure. Points close to the unstable manifold 
of the chaotic saddle have spent a long time in the mix- 
ing region of the flow moving near chaotic orbits with a 
positive Lyapunov exponent. For points precisely at this 
unstable manifold, the backwards trajectories (the ones 
from which the Lyapunov exponent in (|2.12|) should be 
computed) remain in the chaotic saddle, thus leading to 
A > 0. The other trajectories escape from the chaotic set 
in a short time, thus being characterized by a Lyapunov 
exponent equal to zero. 

Thus open flows provide a rather clear example of 
strong space-de pend ence of Lyapunov exponents. Ac- 
cording to Eq. (2.12), the Holder exponent may be dif- 



ferent from 1 only on the unstable manifold of the chaotic 
saddle, thus implying that the transition from smooth to 
filamental structure now only takes place in this fractal 
set of zero measure. The background chemical field is 
always smooth, independently of the value of b. 

To check these ideas, we obtain numerically the dis- 
tributions of chemical fields advected by an open flow. 
Our velocity field is taken from a kinematic model of a 



time-periodic flow behind a cylinder, described in p5[ . 
This flow was found to be qualitatively similar to the 
solution of the Navier-Stokes equation in the range of 
the Reynolds number corresponding to time-periodic vor- 
tex seeding. The concentration pattern shows irregular- 
ities separated by smooth regions (Fig. |[ obtained with 
b = 0.96). This is more clearly observed in the longitudi- 
nal transect. The relation between the singular regions 
and the location of the chaotic saddle can be made patent 
by comparing the gradient of the field with the spatial de- 
pendence of the escape times from the scattering region. 
In particular. Fig. ^ shows the absolute value of the gradi- 
ent of the concentration field |VCoo(i')| that is highly in- 
termittent. It also displays the time (in the time-reversed 
dynamics) that fluid particles initially in a line perpendic- 
ular to the mean flow take to escape the region of chaotic 
motion. Most of the particles leave the region in a short 
time, but longer times appear for initial locations close to 
the chaotic saddle. Clearly, these diverging times are as- 
sociated to the singularities in the gradient distribution. 
By increasing the value of b the flow characteristics (tra- 
jectories, manifolds, escape times, ...) remain unchanged, 
but the singularities in the advected field decrease and fi- 
nally a smooth distribution is obtained. 

A chemical field with the same structure can also be 
obtained in open flows whose time-dependence (and thus 
the chaoticity of advection) is not restricted to a finite do- 
main, by restricting the spatial dependence of the chem- 
ical sources to a finite region. This case was investigated 
in the context of plankton dynamics in |pT| . 

Since the irregularities now appear only on a set of 
measure zero, one could ask if they can have any signif- 
icant effect on measurable quantities. In order to clar- 
ify this, instead of the previous characterization of the 
point-wise strength of the singularities by the Holder ex- 
ponent, let us investigate the scaling of the spatial aver- 
age of the differences SCoo with distance Sr = \Sr\. On a 
one-dimensional transect of unit length the total number 
of segments of length Sr is {Sr)~^ while the number of 
segments containing parts of the unstable manifold (with 
partial fractal dimension D) is N(Sr) Sr~^ . Thus, ac- 



cording to ( ^.12 ) the spatial average of SCcxi along this 
line, ((5Coo(r; Sr)), can be written as 



{SC^{r;Sr)) = {Sr){Sr)-^ (Srf^ + 
+ {Sr)[{Sry^ - {Sr)-^]{Sr) 



(3.1) 



where the first term pertains to the singular, while the 
second one to the smooth component. In the limit (5r — > 
the dominating behavior is 



{SC^{r;Sr)) ^Sr< , 



with 



C = min<j l,l + --i:> 



(3.2) 



(3.3) 
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showing that if D < the average is dominated 
by the smooth component, but if the singularities are 
strong enough (or the fractal dimension of the singu- 
lar set is large enough) they contribute to the scaling of 
((5Coo(r; fr))- The average has been performed along a 
one-dimensional line or transect of the two-dimensional 
pattern. For common velocity fields and transects this 
will be equivalent to the complete average over the whole 
fluid, except in the particular case in which the transect 
is chosen to be completely aligned with the filaments. 



IV. STRUCTURE FUNCTIONS 

The strongly intermittent structure of singularities in 
open flows is an extreme example. There are additional 
inhomogeneities affecting both to the open and to the 
closed flows: although, in the long-time limit the Lya- 
punov exponent is the same for almost all trajectories in 
an ergodic region, deviations can persist on fractal sets 
of measure zero, and as we saw above such sets can con- 
tribute significantly to the global scaling. The origin of 
these inhomogeneities can be traced back by analyzing 
the finite-time distribution of Lyapunov exponents. This 
will be done in the following. For a robust quantitative 
characterization of the filamental structure, accessible to 
measurements, we consider now the scaling properties of 
the structure functions associated with the chemical field. 

The qth order structure function is defined as 



Sg{6r) = (\6Coo{r;Sr)\'^) 



(4.1) 



where ( ) represents averaging over different locations r, 
and 5 is a parameter (we will only consider structure 
functions of positive order (g > 0)). In the absence of 
any characteristic length over a certain range of scales 
the structure functions are expected to exhibit a power- 
law dependence 



Sq{Sr) - 5r'^' 



(4.2) 



characterized by the set of scaling exponents Cq- 

We also note that some of the scaling exponents are 
directly related to other characteristic exponents, such 
as the one characterizing the decay of the Fourier power- 
spectrum r(fc) ^ k~'^ , or the box-counting fractal di- 
mension Dq of the graph of the function Coc{x,y) as a 
function of x by simple relations |l7[ : 



7 = C2 + 1 and D, 



G 



Ci 



(4.3) 



If the Holder expone nt of the field has the same value 
everywhere, given by (2.12), the scaling exponents of the 
resulting mono-ajfine field are simply 



(4.4) 



(we have assumed b < X). In the case of our chemical field 
this, however, would only hold in the non-generic case 



in which the stretching by the flow is spatially uniform 
(e.g. like in a simple area-preserving baker's map). In 
general, the singular spatial inhomogeneities of the Lya- 
punov exponent could be understood by realizing that 
the finite-time stretching rates, or finite-time Lyapunov 
exponents fl^ , have a certain distribution around the 
most probable value. This distribution approaches the 
time-asymptotic form fl^ , [l7|| : 



(4.5) 



where G{X) is a function characteristic to the system 
(velocity field, in this case), with the property that 
G{Xq) = G"(Ao) = and G(A) > 0, where Xq is the 
most probable v alue of the infinite-time Lyapunov expo- 
nent. The form (4^ ) is valid only for hyperbolic systems. 
Non-hyperbolicity can strongly affect the distribution at 
small values of A but around Aq and for larger values 
it remains a good approximation. As we shall see later 
only this region contributes to the structure functions of 
positive order. 

As time increases the distribution becomes more and 
more peaked around Aq. The (Lebesgue) area of the set 
of initial conditions with finite-time Lyapunov exponents 
in a small interval (A, A + SX) that excludes Aq decreases 
at long times with a dominant exponential behavior: 



SAx{t) - e~^(^)*(5A 



(4.6) 



showing that only sets of measure zero can have Lya- 
punov exponent different from Aq in the t —^ 00 limit. 
Such sets, however, can still have nonzero fractal di- 
mensions. At finite times, the area (4.6) encloses the 
final anomalous set, with a transverse thickness that, 
due to stretching by the chaotic advection, decreases like 



hit) 



-At 



The number of boxes needed to cover the 



set of area SA\ (t) using boxes of size l\ (t) is 



Nxit) 



g[2A-G(A)]t _ ^(^)i^ 



(4.7) 



that gives the dimension for the set to which this area 
converges in the infinite-time limit: 



D{X) = 2- 



GjX) 
X ' 



(4.8) 



Thus, an arbitrary line across the system will be found 
composed by subsets of dimension D{X) = D{X) — 1 each 
one characterized by different values A of the Lyapunov 
exponents and in consequence of the Holder exponents 
a{X) = min{6/A, 1}. 



Now, the scaling exponents in (L2) can be readily ob- 
tained. The number of segments of size 6r belonging to 
a subset characterized by Lyapunov exponent A scales as 
A^(A) ~ 5r^^'^^\ while the total number of such non- 
overlapping segments scales as ~ Sr^^. Thus, the struc- 
ture function can be written as 
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5r2-^W|^Coo(r(A),fr)|''dA 



Sr^-DiX)g^,^^^ T"" ^^2-D(A)^^,6Arf;^ (4 9) 
A,„i„ Jb 



In the limit (5r — s- the integrals are dominated by a 
saddle point and, after some manipulations, the scaling 
exponents in ( |4.2D are obtained as 

(4.10) 

The right hand side can be seen as a family of lines in 
the q — Cq plane labeled by the parameter A so that the 
value of Cg is given by the lower envelope of these lines. 
Note that the shape of G(A) for A small enough becomes 
irrelevant for determining because of the minimum 
condition. Thus multifractality, characterized by nonlin- 
earity in the g-depencence of Cg, is affected onl y by the 
largest stretching r ates in the flow. Equation (3^) is a 
particular case of ( 4.10| ) for q = 1 and in the approxi- 
mation of considering a single value of A on the chaotic 
saddle. 

According to (|4.1C|) the gth order structure function 



is dominated by a subset characterized by a Lyapunov 
exponent Xq. Applying the extremum condition to (4.1C) 
we obtain an equation for Xq 



GiXq) 



(4.11) 



that can be substituted into ( 1.10| ) to obtain the qth order 
scaling exponent. 

We have analyzed numerically the c hemical decay un- 
der advection by the closed flow ( ^.13 ) to check the the- 
oretical predictions above. Numerically computed his- 
tograms of the finite-time Lyapunov exponents are shown 
in Fig. ^ and the corresponding G(A) functions are repre- 
sented in Fig. p. The G(A) functions obtained from his- 
tograms corresponding to different times coalesce except 
in the small-A region, thus confirming that ( [4.5D correctly 
describes the observed distribution. 

Numerically calculated scaling exponents (i.e. ob- 
tained by direct application of Eq. (4.1) averaging over 
many one-dimensi onal t ransects) , and the family of lines 
corresponding to ( 4.10| ) based on the histogram of the 
finite-time Lyapunov exponents of Fig. |^ are shown in 
Fig ^, where the prediction of the mono-fractal approx- 
imation (Cg = qb/Xo) is also shown. The mono- fractal 
approximation appears to be accurate for small q. The 
graph-fractal dimension or the widely used Fourier power 
spectrum exponent are related to Ci and C2 by Eqs. ( |4.3| ) 
so that their estimate based on the mono-fractal descrip- 
tion that considers just the bulk value of the Holder ex- 
ponent can deviate from the actual values. 

In a recent work by Nam et al. the power spectrum 
of a decaying scalar field (with space-dependent decay 



rate) has been investigated and related to the distribu- 
tion of finite-time Lyapunov exponents of the advecting 
fiow. The result for the spectral exponent obtained in 
|l8| using an eikonal-type wave packet model |jl^, and 
taking in to ac count finite diffusion, is consi sten t with our 
formula ( 4.10 ) (that for g = 2, and with (L3) gives the 
value of the spectral slope) obtained in the non-diffusive 
limit. 

The function G(A) is characteristic to the advecting 
fiow. Let us now consider a special case where we ap- 
proximate G(A) by a parabola 



G(A) 



(A-Ao)^ 
2A 



(4.12) 



This can be thought as the first term in a Taylor expan- 
sion around Aq , which is a good approximat ion to obtain 
the small-g scaling exponents. In this case ( 4.11 ) can be 
solved explicitly 



Ag = V(Ao)2 + 2(z&A 
This gives the scaling exponents 



'Ao 



2qh Ao 



(4.13) 



(4.14) 



The above relation has been obtained recently by 
Chertkov in , where the problem of advection of de- 
caying substances was considered in a probabilistic set- 
up, using stochastic chemical sources and a random ve- 
locity field that is spatially smooth but delta-correlated 
in time. The distributi on of stretching rates was assumed 
to be Gaussian as in (4.12). This assumption could be 
realistic in many cases, and could give good estimates 
for the scaling exponents for small q. For higher-order 
moments, however, higher-order terms in the expansion 
of G(A) can become important. Moreover, the possible 
values of A could be limited by a finite maximum value 
AmaKj e.g. in time-periodic fiows, where the finite time- 
Lyapunov exponents cannot have arbitrarily large values. 
This implies that the scaling exponents for q > q* , where 
A, 



9* 



Xmax, should display a simple linear dependence 



qb + G{Xmax) 
^max 



(4.15) 



that differs from the q-i behavior of (4.14) for large q 



V. SUMMARY AND DISCUSSION 

Multifractality of advected fields generated by chaotic 
advection has been observed previously in the case of 
passive advection with no chemical activity {h — 0) pl| . 
It was shown that the measure defined by the gradients 
of the advected scalar field has multifractal properties 
and its spectra of dimensions Dq has been related to 
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the distribution of finite-time Lyapunov exponents 
This multifractahty, however, does not affect the slope of 
the power spectrum [^,^ (the so cahed Batchelor spec- 
trum: r(fc) ^ k~^) or scahng exponents of the s truct ure 
functions (Cg = for all q, as it can be seen from (4.1C) by 
setting b = Q). The effect of multifractahty on the power 
spectrum has a character transient in time, moving to- 
wards smaller and smaller scales and finally disappearing 
when reaching the diffusive end of the spectrum. In the 
stationary state only the diffusive cut-off of the power 
spectrum is affected |24| that can still be important for 
the interpretation of some experimental results. By com- 
paring these results for the conserved case with the ones 
presented here for the decaying scalar we can conclude 
that, although the origin of the multifractahty is the same 
in both situations -the non-uniformity of the finite-time 
Lyapunov exponents- in the presence of chemical activ- 
ity this has stronger consequences (non-Batchelor power 
spectra and anomalous scaling). 

We have presented a simple mechanism that can gener- 
ate multifractal (or more precisely, multi-affine) distribu- 
tions of advected chemical fields. The main ingredients 
are chaotic advection and linear decay of the advected 
quantity in the presence of non-homogeneous sources. 
Essentially the same mechanisms were considered in [|oj , 
but with stochastic time-dependencies both in the flow as 
in the chemical sources, considering advection by the spa- 
tially smooth limit of the so-called Kraichnan model gen- 
erally intended to represent turbulent flows. Our results 
stress that anomalous scaling may appear in simple regu- 
lar (e.g. time-periodic) laminar flows where stochasticity 
appears just as a consequence of the low-dimensional de- 
terministic chaos generated by the Lagrangian advection 
dynamics. In addition we have provided numerical evi- 
dence for the theoretical predictions. 

Chaotic advection is characteristic to most time- 
dependent flows. The linear decay of the advected sub- 
stance is in fact just the simplest prototype of a family 
of chemical reaction schemes, where the local dynam- 
ics converges towards a fixed point of the chemical rate 
equations. The local dynamics can also be generated by 
non-chemical processes, e.g. by biological population dy- 
namics in the case of plankton advection , or by the 
relaxation of the sea-surface temperature towards the lo- 
cal atmospheric value |0 . Inhomogeneities of the chem- 
ical sources or of other parameters of the local dynam- 
ics arise naturally in these contexts so that we expect 
our results to be of relevance in biological and geophys- 
ical settings. In fact, fractality and multifractahty have 
been already observed in these contexts, for example in 
the distribution of stratospheric chemicals (e.g. ozone) 
p5| , p6| , and in sea-surface temperature and phytoplank- 
ton populations The structure of these fields has 

been sometimes associated with turbulence of the advect- 
ing flow. We think that the simple mechanism, able to 
generate complex multifractal distributions, investigated 
in this Paper can be at the origin of some of the struc- 
tures observed in geophysical flows. Further work in this 



direction could help on the interpretation of geophysical 
data in order to gain quantitative information about the 
processes involved. Laboratory experiments seem also to 
be feasible. 



ACKNOWLEDGMENTS 

Helpful discussions with Peter Haynes and Oreste Piro 
are greatly acknowledged. This work was supported 
by CICYT (Spain, project MAR98-0840) and DGICYT 
(project PB94-1167). Z.N. was supported by an Euro- 
pean Science Foundation/TAO [Transport Processes in 
the Atmosphere and the Oceans) Exchange Grant. 



[1 
[2: 

[4: 
[s: 
[6: 

[7 

[8 
[9 
[10 

[11 

[12 

[13 

[14 

[15 
[16 

[17] 



J.M. Ottino, The Kinematics of Mixing: Stretching, 
Chaos, and Transport (Cambridge Univ. Press, Cam- 
bridge, 1989). 

H. Aref, J. Fluid Mech. 143, 1 (1984); J.H.E. Cartwright, 
M. Feingold, and O. Piro in Mixing: Chaos and Turbu- 
lence, H. Chate and E. Villermaux Eds. (Kluwer, Dor- 
dretch, 1999). 

I. R. Epstein, Nature 374, 321 (1995). 

S. Edouard, B. Legras, F. Lefevre and R. Eymard, Na- 
ture, 384, 444 (1996). 

C. Marguerit, D. Schertzer, F. Schmmitt, S. Lovejoy, J. 
Mar. Sys. 16, 69 (1998). 

G. Metcalfe and J.M. Ottino, Phys. Rev. Lett. 72, 2875 
(1994). 

Z. Toroczkai, Gy Karolyi, A. Pentek, T. Tel and C. Gre- 
bogi, Phys. Rev. Lett. 80, 500 (1998); Gy. Karolyi, A. 
Pentek, Z. Toroczkai, T. Tel, and C. Grebogi, Phys. Rev. 
E 59, 5468 (1999). 

Z. Neufeld, C. L opez and P.H. Hayn es, Phys. Rev. Lett. 
82, 2606 (1999). [chro-dyn/9906019 

J.-P. Eckmann and D. Ruelle, Rev. Mod. Phys. 57, 617 
(1985). 

E. Abraham, A note on the wavenumber spectra of sea- 
surface temperature and chlorophyll, preprint (1999). 
C. Lopez, Z. Neufeld, E. Hernandez-Garci'a, P.H. Haynes, 
Chaotic advection of reacting subst ances: Plankton dy - 
namics on a meandering jet, preprint chao-dyn/9906029. 



E. Ott, Chaos in Dynamical Systems (Cambridge Univ. 
Press, Cambridge, 1993). 

T. Tel, in Directions in Chaos, Vol. 3, edited by Hao 

Bai-Lin (World Scientific, Singapore, 1990). 

E. Ziemniak, C. Jung and T. Tel, Physica D 76, 123 

(1994) 

C. Jung, T. Tel and E. Ziemniak, CHAOS 3, 555 (1993). 
J.C. Sommerer, H.-C. Ku and H.E. Gilreath, Phys. Rev. 
Lett. 77, 5055 (1996). 

T. Bohr, M. Jensen, G. Paladin, A. Vulpiani, Dynami- 
cal Systems Approach to Turbulence (Cambridge Univ. 
Press, Cambridge, 1998). 



8 



[18] 
[19] 
[20] 
[21] 
[22] 
[23] 

[24] 
[25] 
[26] 
[27] 



K. Nam et al., The k spectrum of Bnite lifetime passive 
scalar in Lagrangian chaotic fluid flows, (preprint, 1999). 
T. Antonsen, Jr. Z. F. Fan and E. Ott, Phys. Rev. Lett., 
75, 1751 (1995). 

M. Chertkov, Phys. Fluids 10, 3017 (1998). 



chao-dyn/9803007 



E. Ott, T. M. Antonsen, Jr., Phys. Rev. Lett. 61, 2839 
(1988). 

T. M. Antonsen, Jr. and E. Ott, Phys. Rev. A, 44, 851 
(1991) 

M. Falcioni, G. Paladin and A. Vulpiani, Phys. Rev. 
Lett., 64, 698 (1990); E. Ott, T.M. Antonsen Jr., Phys. 
Rev. Lett., 64, 699 (1990). 

G.-C. Yuan et al.. Power spectrum of passive scalars in 

two-dimensional chaotic flows (preprint, 1999). 

J. Bacmeister et al, in Gravity wave processes, edited by 

K. Hamilton (Springer- Verlag, New York, 1997). 

A.F. Tuck and S.J. Hovde, Geophys. Res. Lett. (1999, to 

appear). 

L. Seuront, F. Schmitt, Y. Lagadeuc, D. Schertzer, S. 
Lovejoy and S. Frontier, Geophys. Res. Lett. 23, 3591 
(1996). 



0.0 0.2 0.4 0.6 0.8 1.0 



0.05 - U,!" % 

Ha 



-0.05 
-0.10 



0.0 0.2 0.4 0.6 0.8 1.0 




FIG. 2. Top: a snapshot of the chemical concentration Coo 
for 6 = 0.1 and A — 2.67, so that a filamental structure is 
obtained. The lower panel shows a horizontal cut along the 
line y = 0.25, clearly displaying the fractal nature of the field. 
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FIG. 1. Top: a snapshot of the chemical concentration Coo 
obtained for b — 4.0 and A = 2.67, so that a smooth distri- 
bution is obtained. Bottom: a horizontal cut along the line 
y = 0.25. 



FIG. 3. Top: a snapshot of the chemical pattern Coo 
formed in the wake of a cylinder (the black semicircle at the 
left is half of its section). Mean fiow is from left to right. We 
have used the streamfunction given in reference [15]. b = 0.96, 
and all the streamfunction parameter values are the same as 
in [15] except the boundary-layer thickness of the cylinder 
which here takes the value a = 20.0, and the vorticity strength 
which is Lj = 35.06 in our calculations. Bottom: a horizontal 
cut taken along the line y = 1.0. 
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FIG. 4. Top: The absolute value of the gradient of the 
chemical field in Fig. ^ along the line x = 7.3. The lower figure 
shows the escape time for particles along this line (x — 7.3). 
This is calculated by computing the time that every single 
particle takes to arrive to the line x = —2.0 (far from the 
chaotic wake region) in the backwards-in-time dynamics. 
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FIG. 6. The fun ction G(A), obtained from the distributions 
in Fig. ^ and Eq. (4.5). Collapse of data for the three times 
into the same curve confirms the validity of ( [l.5| ) . The dotted 
line is a parabolic fit (0.4(A — 2.67)^) that provides a good 
approximation near the minimum. 
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FIG. 5. The distribution of finite-time Lyapunov expo nent s 
at three different times, obtained for the closed fiow (2.13) 
with T = 1.0 and U — 1.2. The long-time mean Lyapunov 
exponent is Aq — 2.67. 
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FIG. 7. The scaling exponents (^q for the situation of Figs. ^ 
and ^ and with b — 1.0. Thick line: the mono-fractal ap- 
proximation (q = qb/Xo. Thin lines: the curves (q = q and 
Cq ~ [<lb + G{X)]/ X, for different values of A; the numeric al va. 1- 
ues of G(A) are obtained from Fig. H According to Eq. ( [l.lO| ), 
the actual values of the scaling exponents are given by the 
lower envelope of this set of curves. This is confirmed by the 
numerically determ ined values of (^q (crosses). Dashed line: 
the approximation (4.14). 
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